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ABSTRACT 



H \ A previous analysis of two-body Cabibbo allowed nonleptonic decays of D^ mesons and of 

Cabibbo allowed and first-forbidden decays of D^ and D~^ has been adjourned using more 
recent experimental data and extended to the Cabibbo forbidden decays of D^ . Annihila- 
tion and W-exchange contributions as well as final state interaction effects (assumed to be 
dominated by nearby resonances) have been included and are in fact crucial to obtain a rea- 
sonable agreement with the experimental data, which show large flavour SU(3) violations. 
New fitting parameters are necessary to describe rescattering effects for Cabibbo forbidden 
D^ decays, given the lack of experimental informations on isoscalar resonances. We keep 
their number to a minimum - three - using phenomenologically based considerations. We 
also discuss CP violating asymmetries. 



'■*'' partially supported by the European Community under the Human Capital and Mo- 
bility Programme, contract CHRX-CT93-0132. 



This paper replaces the previous one with the same title submitted to the ICHEP94 con- 
ference (ref. gls0658), also circulated as ROMA preprint n.1026-1994 / NAPOLI DSF-T- 
18/94. The 1994 edition of the Particle Data (Phys. Rev. D 50 (August 1994) part I) 
did in fact show noticeable changes in the experimental branching ratios for D decays. We 
therefore made a new analysis, that we report here. 

1. Introduction 

A theoretical description of exclusive nonleptonic decays of charmed hadrons based 
on general principles is not yet possible. Even if the short distance effects due to hard gluon 
exchange can be resummed and an effective hamiltonian has been constructed (recently, at 
next-to-leading order |^ ) , the evaluation of its matrix elements requires nonperturbative 
techniques. Waiting for future progress in lattice QCD calculations one has to rely on 
approximate methods and/or models. 

The largely different lifetimes of charmed hadrons make it clear that the infinitely 
heavy quark limit is quite far from the actual situation. Therefore, the expansion in in- 
verse powers of the heavy quark mass characteristic of heavy quark effective theory 
is presumably not a useful tool in this case. Moreover, the methods of HQET are not 
obviously extended to cope with exclusive hadronic decays. On the other hand, the simple 
factorized ansatz for the matrix elements is known not to describe properly Cabibbo al- 
lowed D^ decays. The color-suppression of some contributions seems in fact to be stronger 
than the factor 1/3 expected from QCD and the data exhibit large phase differences 
between amplitudes with definite isospin. We are thus forced, still using the factorization 
approximation as a starting point of the matrix element evaluation, to include important 
corrections due to rescattering effects in the final states. This we do assuming the domi- 
nance of nearby resonances and taking from experiment - when possible - their masses and 
widths. We also include W-exchange and annihilation contributions that turn out to be 
larger than generally believed. The presence of nearby resonances may well have the effect 
of increasing these terms relative to their naive PCAC estimates. 

In two previous papers the Cabibbo allowed [^ two-body decays of charmed mesons 
were described in the framework discussed above and the model was applied to the analysis 
of Cabibbo first-forbidden decays of the charged mesons ||] D"*" and Df and to their CP 
violating asymmetries. 

The considerable success of that analysis prompts us to extend it to the Cabibbo 
forbidden two-body decays of D^. The recent experimental determination of the branching 
ratio B.R.(D+ -^ Tt+n^) = 0.25±0.07 % 0] - that agrees with our prediction 0] - allows 



to perform an amplitude analysis on the complex of D ^ nn decays, that shows a large 
(~ 90°) phase difference between / = and 1 = 2 amplitudes. Moreover, a comparison 
of the 1 = 2 amplitude with the I = ^ from D~^ — > ti~^K shows a considerable violation 
of flavour SU(3) in the direction of larger tttt amplitudes; on the other hand, it is known 
since a long time and recently confirmed that the ratio of D^ decay branching fractions 
to K~^K~ and to tt+tt" is much larger than the SU(3) prediction (i.e. 1), showing an 
opposite pattern of SU(3) breaking in exotic and nonexotic channels. Another striking 
signal of the importance of SU(3) violations is given by the value, quite similar to other 
Cabibbo forbidden decays, of the B.K.{D^ — > K K^) that should be vanishing in the 
symmetric limit. 

Our model describes satisfactorily the experimental situation. For what concerns 
SU(3) breaking in exotic channels it is the combination of several small effects that yields 
a large result. These effects would not be enough to explain the large K^K~ to tt+tt" 
ratio and to produce a nonvanishing decay rate to K K^: in the nonexotic channels the 
rescattering effects are essential. 

In order to introduce these rescattering effects we need to know masses, widths and 
couplings of yet unobserved spinless, isoscalar resonances with positive and negative parity 
and masses around 1.9 GeV. One expects for each parity two resonances of this type, a 
SU(3) singlet and a member of an octet, that generally mix among themselves. Such a 
large number of new parameters to fit eight new data (or limits) for branching ratios is 
obviously unappealing, unless some arguments can be given to reduce it. In the following 
we will show that reasonable phenomenological assumptions may reduce the number of 
new parameters to three. 

We have to determine these by a fit to the data. Before doing that we repeated the 
fit to all Cabibbo allowed and to charged meson first-forbidden decay branching ratios, 
which in the meantime have got lower error bars and in some cases have also changed. 
The model is therefore passing a more demanding test. 



2. Decay amplitudes in the factorized approximation 

The effective weak hamiltonian for Cabibbo allowed nonleptonic decays of charmed 
particles is given by {Uij are elements of the CKM matrix) 

H^fr""' = % U-d K, [C^ S^ 7/. (1 - 75) Ca U^ 1^ (1 - 75) rf/3+ ^ ^ 

V2 (2.1) 

Ci w" 7^ (1 - 75) Ca s^ -i^ (1 - 75) dfi\ + h.c. , 
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while for AC = —AS" processes the hamiltonian is obtained from the same equation 
with the substitution s ^->- d. The effective weak hamiltonian for Cabibbo first-forbidden 
nonleptonic decays is 

C 

+ -Jk ^«s ^cs [Ci Qi + C2 Q2] - (^2.2) 



V2 
V2 



- ^UubU:^ X^ C,Q, + h.c. 



i=3 



-^a: 



In eq. ( p7^ ) the operators are [||: 

Qi = u^ 7^ (1 - 75) dfi d^ r (1 - 75) cc 

Qi = ^" 7^ (1 - 75) d^ d" Y (1 - 75) c^, 

Q3 = ^^" 7a. (1- 75) Ca Y. ^^7^(l-75)g/?, 

<? 

^4 = ^" 7,. (1 - 75) cp J2 ^^ ^'^ (1 - 75) ga, (2.3) 

Q5 = ^i"7M(l-75)ca 5^ ?^7^(l+75)g/3, 
Q6 = ^^"7a.(1-75)c/3 ^ ?''7^(l+75)gc.. 

The operator Qf (QI) in eq. ( p.2|) is obtained from Qf ((52) with the substitution 
((i ^ s). In eqs. (|2.1J ) and (|2.3| ) a and /3 are colour indices (that we will omit in the 
following formulae) and in the "penguin" operators q [q) is to be summed over all active 
fiavours (m, d, s). 

If we neglect mixing with the third generation {Uub = and Uus U*g = —Uud U*^ = 
sinOc cos 9c) then the three effective hamiltonians 

ttAC=AS ttAS=0 ttAC=-AS 

cos2^c ' 72 sin ^c COS ^c ' sin^ ^c 

form a C/-spin triplet. Therefore, in the limit of exact fiavour SU(3) symmetry a number of 
relations between decay amplitudes should hold. We shall discuss some of them in Section 
4 and we will see that they are often violated rather strongly. 

We have evaluated the coefficients Ci at the scale 1.5 GeV using the two- loop anoma- 
lous dimension matrices recently calculated by A. Buras and collaborators [|l|, assuming 



\MS _ 3QQ MgY This value, which corresponds to the best agreement between the ex- 
perimental data and the theoretical results on the exclusive decay channels of D mesons, 
is compatible with the experimental determination from LEP measurements [§ . The coef- 
ficients at next-to-leading order are renormalization scheme dependent: we assume in the 
following the values obtained using the "scheme independent prescription" of Buras et al., 
namely Ci = -0.628, C2 = 1.347, C3 = 0.027, C4 = -0.057, C5 = 0.015, Cq = -0.070. 

In the factorized approximation the matrix elements of -ffe// are written in terms 
of matrix elements of currents, {V5)^ = q' 'j^ q and (A^,)^ = q' 'j'^'J^q. 

We recall the definitions of the decay constants for pseudoscalar (tt, K. . .) and 
vector (p, K*. . . ) mesons 



<P,(p)|AnO>=-i/p,p^ , 
<V,{p,X)\V^\0>=MUv.e*^{X), 

and the usual definitions for the matrix elements of the currents: 

Mf -Mf 

r 



(2.5) 



(2.6) 



f * . a M^ — M^ 

< Vi\V^\P, >=2V{q') Z ^ ^;/ . 

' ' ^ ^ ^ Mj + Mi 

To avoid the presence of a spurious singularity at q'^ = 0, one has to require that 

/+(0) = /o(0)^^,,„ (2.7) 

, , , Mi + M. , , Mi- Mi ^ , , , , 



The semileptonic decay rate for D^ -^ K e'^v ilO] indicates a value Vcs — 0.79; 



assuming SU(3) symmetry for the weak charges, we will set 

Vcs = Vcd = Vcu = 0.79 (2.9) 
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Using the definition (|2.8|) and the data on the semileptonic decay D^ -^ K*~e~^u [^ 
we obtain acs — -54 ± .13 (E653 CoUaboration) or Ucs — -71 ± .16 (E691 CoUaboration) . 
Different lattice QCD calculations |jl^ give similar results: in average acs — -74 ± 0.15. 
The limited statistics for the data on D^ -^ p~e^u decay does not allow an analysis of the 
different form factors, but within large errors the measured branching fraction is larger than 
theoretical predictions based on quark models |]T^. Lattice results and results obtained 
in the framework of QCD sum rules by using ( |2.8D on SU(3) breaking are inconclusive 
because of large errors on the A2 form factor. 

Since the data on D meson decays show large SU(3) breaking effects and since the 
axial charges are not protected by the so-called AdemoUo-Gatto theorem, in our fit we 
allowed acs and acd = cicu to vary between .5 and 1 independently. The values chosen by 
the fit are acs = 0.59, consistent with experimental data, and acd = 1 (in fact, an even 
better fit would be obtained allowing larger values for acd)- We note that these values do 
not agree with the direct QCD sum rule calculation oi AQ[q^) performed in |]I^, where the 
authors conclude that SU(3) breaking effects should be small {acs/cicu = 1-10 ± 0.05) and 
the q^ dependence of the form factors compatible with a polar dependence dominated by 
the 0~ pole. 

The q^ dependence of the form factors is assumed to be dominated by the nearest 
singularity. This entails for the c-quark decay terms the usual simple-pole form factors 



„ ,2\ ^_ '^cu(s,d) 



Mq 



2^ 



(2.10) 

Vcuis,d) ^ ' 



1 - q'iMl, 



'^(.)(0+) 



and analogous expressions, with the mass of the lightest particle with appropriate quantum 
numbers, for the other form factors. 

In the W-exchange and annihilation terms, however, the large and time-like (^ 
values needed, together with the suggested existence of resonances with masses near to 
the D-meson mass, make a prediction based on the lightest mass singularity unjustified. 
These terms depend on the matrix elements of current divergences between the vacuum 
and two-meson states. We write them, with the help of the equations of motion, in the 
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(2.11) 



way indicated in the following examples: 

< K-n+\d''{Vf)^\0 >= i {rus - ma) < K'n+l sd |0 >= 

JVf2 

Id 

< K-p+\d''{Ai)f,\0 >= i (m, + rud) < K'p+l s-f^d |0 >= 

Id 



We assume SU(3) symmetry for the matrix elements of scalar and pseudoscalar densities, 
and express all of them in terms of Wpp, Wpv- In our approach the VF^'s are free param- 
eters of the fit. Their magnitude turns out to be considerably larger than what one would 
obtain assuming form factors dominated by the pole of the lightest scalar or pseudoscalar 
meson, i.e. K^{U30) or K{497). 

We note that to obtain the amplitudes for Cabibbo first-forbidden decays one has to 
evaluate matrix elements like < 77|J7^75(i|0 > or, for penguin operators, < r]\dj5d\0 >. To 
get the correct result it is necessary to take into account the anomaly of the singlet axial 



current; we have followed the method discussed in [15|. In this scheme the 77 — 77' mixing 
angle, 9nri', results to be equal to —10°. Remarkably, this value which is consistent with 
the Gell-Mann Okubo mass formula, is also perfectly compatible with the experimental 



value of r(77 -^ 77) obtained by two-photon production experiments [16|. Therefore the 



rj ~ Tj' mixing angle is not a parameter of the fit, as it was in our previous analyses. 

If the final K meson in Cabibbo allowed decays is neutral, it has been observed as 
a short-lived neutral K, Kg- There is therefore an interference between Cabibbo allowed 
{D ^>- K + X) and doubly suppressed {D -^ K^ + 
the experimental data Texp [0 using the definition 



{D ^ K + X) and doubly suppressed {D -^ K^ + X) decay amplitudes. We have fitted 



V{D ^Ks + X) = ]^ re.p(D ^ kV X) , (2.12) 



and included FSI modifications also for the doubly Cabibbo suppressed part of the ampli- 
tude. The correction due to this effect is not negligible and it helps in obtaining a better 
fit to the experimental data. 

We write now a few examples of amplitudes for Cabibbo allowed and first forbidden 



decays in the factorized approximation: 

A^{D^ -^ K-n+) = 
Gp 



V2 



+ {Ci+^C2) <K-n+\{Vf)^\0><0\{Air\D'>>]= (2.13) 



^^UudU:, [iC2 + ^Ci)U <K-\d^'iVn^\D^> + 



V2 

+ iC,+CC2)fD <K-n+\d^iVf)^\0>]. 

In ( |2.13|) the two terms correspond to c-quark decay and W-exchange, respectively . 

AUD+ ^ Ksp+) = 

+ {Ci+^C2) < lf\{Ai),\0 > < p+\{Air\D+ >] + 

+ {C^+iC2) <K'^\{A2^Q><p+\{AlY\D+>]= (2.14) 
= -^U^dU:, [{C2 + ^C,)f,M,e*^{p+) <lf\{V:),\D+> + 
+ {Ci+^C2)fK <p+\d^{Al)^\D+>] + 
+ ^U^sU:, [{C2+^C,)fo <K'p+\d^{A:),\0> + 
+ {Ci+^C2)fK <p+\d^{Al)^\D+>]. 

In (|2.14| ) we have given an example of an amplitude with K and K^ interfering contribu- 
tions, that also contains an annihilation term. 

A^{D+ -^ 7r%+) = 

= -^U^dU:, [{02+^0 <AiVdl,\D+><7r+\{Atr\0> + 

+ [Ci +^C2) < 7:+\{V:),\D+ > < n\Air\0>] = 

-{c, + iC2)^ <7r+|a^(Kr),|D+>] = 

— I ^ TT .TT* (n. I r<^\ (^ __l c\ f ^ ^+i;^m/'t/c^ i n+ 



u^dKd (^1 + ^2) (1 + 0/^ < 7r+|a^(K^) ji?+ > . 



A^iD'^ -^ 7r-7r+) = 
= -^U^^U:, (C2 + e^i) <n-\{yi)^\D'^><n+\{AiY\Q> + 

+ ^U^,U:, [{C^ + ^Cs) <n-\{Vi)^\D^><n+\{Atr\0> + 

-2(C6 + ^C5) < 7r-7r+|™|0 > < 0|w75c|L'° > + 
+ 2(C6 + ^C5) <7t-\dc\D'^ > <n+\u-f5d\0>] = 



■72 

Gf 



(2.16) 



+ 2z(C6 + eC5) <7r-7r+|m|0> -^ ^ + 

ruu + rnc 

-2z(C6 + eC5) <iT-\dc\D'^> -^ — ^1. 

In (|2.15| ) and ( |2.16| ) we give examples of Cabibbo forbidden decay amplitudes, to the second 
of which penguin operators contribute. 

The parameter ^ appearing in the above equations should be equal to | in QCD. 
Since however other colour suppressed, nonfactorizable contributions have been neglected 
we will consider it as a free parameter to be fitted to the experimental data, following ||^. 
The result of the fit favours a value ^ — 0. 



3. Final state interaction effects 

We make the assumption that FSI are dominated by resonant contributions, and 
we neglect the phase-shifts in exotic channels. In the mass region of pseudoscalar charmed 
particles there is evidence, albeit not very strong [Q, for a J^ = 0~ K{1830) (with F = 250 
MeV and an observed decay to K(p ^) and a J^ = 0" 7r(1770) with F = 310 MeV |l|]. 
The coupling of an octet of 0~ P resonance to 0~ + 1~ {PV) channels is determined from 
charge conjugation and SU(3) symmetry to be: 



hfabc{d^Pa)V,^Pc. (3.1) 

Equation (|3.1| ) implies: 

nK ^ PV) _ T..AK) .32. 

r(^ - PV) - "-^ - T^^,in) ' ^'^■'> 

consistent with the assumption that the P resonances decay predominantly into the lowest 
lying P and V mesons, which we shall make for simplicity. 
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For Cabibbo allowed and doubly forbidden D^ decays and for Cabibbo first and 
doubly forbidden D^ decays, the FSI effect modifies the amplitudes in the following way: 

A{D -^ Vh Pk) =A^{D^ Vh Pk) + c;,fc[exp(z58) - 1] $^ Ch,k, A^{D -^ Vk, Pk,) , (3.3) 

h/kf 

In ( |3.3| ) Chk are the normalized {J^^hk ~ 1) couplings PPV and 

T(P) 

sinSs expfz^s) = ^^. (3.4) 

2{M~-Md)-iT{P) 

where P is the resonance appropriate to the decay channel considered (tt or K) . Equation 
( p^ ) determines Ss up to a 180° ambiguity. The choice can be made according to the 
number of resonances and bound states [^ that are present in the channel at lower 
energies, each one of them increasing the phase- shift by tt. In this case the P resonance 
may be assumed as the third resonance in the channel, with 5s close to ^vr (in our recent 
preprint a different choice was made for the PV channels). 

For Cabibbo forbidden L>° —>■ PV decays one expects rj and ry' resonances to take 
also part to rescattering. Due to charge conjugation invariance, the singlet components 
have vanishing coupling and the combined effect of the two expected isoscalar resonances 
may be described by a phase attached to the isosinglet octet part of the decay amplitudes. 
This phase is the only added parameter for these channels to be fitted. The fitted value is 
~ 243°, corresponding to two resonances with masses one below and the other above the 
D^ mass. 

Coming now to FSI effects for parity violating D —>■ PP decays, we note that 
some evidence exists for a J^ = 0"*" resonance Kq (with mass 1945ibl0ib20 MeV, width 
201±34±79 MeV and 52±14% branching ratio in Kir ^^). No ao isovector resonance has 
been observed up to now in the interesting mass region. In |^ we assumed its existence 
and we estimated its mass from the equispacing formula 

Ml=M],,-Ml + Ml (3.5) 

In the fit we allowed the mass and the width of Kq to vary within the experimental bounds. 



The best fit values are 1928 MeV and 300 MeV, respectively. From Eq. (3^) we get M^q = 
1869 MeV. 

The SU(3) and C invariant coupling of a scalar S resonance to two pseudoscalar 
mesons is: 

-^ 9888 dabc Pa Pb Sc+g818 { Pa Pq + Pq Pa) Sa+ , . 

+fi'881 Pa Pa So+9lll Pq Pq Sq- 

10 



In ( |3.6|) a,b,c= 1 ... 8 are SU(3) indices, Pq and 5*0 are SU(3) singlets. 

Assuming that the S resonances decay dominantly to a pair of mesons belonging 
to the lowest mass nonet, one obtains from (|3.6| ) the branching ratio for (Kq — > Kn) as 
a function of the ratio of coupling constants r = gsis / 9888- The further assumption of 
nonet symmetry would imply r = 1. The experimental data allow two possible values for 
r, one positive (and consistent with 1) and another negative and close to —1. The best fit 
value of r is —0.84, corresponding to a branching ratio for the decay of the Kq resonance 
in Ktt of about 64%. 

The description of rescattering effects for Cabibbo forbidden D^ decays is compli- 
cated by the presence of yet unobserved /o and /q isoscalar resonances, which should be 
singlet-octet mixtures. Denoting by |/o) the lower mass state, we define 

l/o) = sine/) l/s) +cos(/) l/i) 
|/^)=-cos(/>|/8)+sin</>|/i). 



The results will also depend on the parameters a = gs8i / 9888 and c = gm / 9888, see ( |3.6|) . 
To reduce the number of new parameters, we assume that these scalar resonances 
behave similarly to the tensor mesons /2 (1270) and /2 (1525): the /2 is very weakly 
coupled to TTTT, and the /2 has in turn a small coupling to KK. In order to forbid the 
/q -^ TTTT decay, we required a and the mixing angle (f) to be related by 

1 



a/2 tani 



The value tancj) = v2 would then imply a vanishing branching ratio for the decay /o -^ 
KK. The best fit value is tanc/) = 1.14, not very far from 1.41. 

The fit has therefore two new free parameters, A^ = M?, — M? and (p. For any 
pair (A^, (p) there are two possible values for c, that are solutions of a quadratic condition 
coming from the requirement of unitarity of the rescattering transformation. 



4. Comparison with experimental data on Branching Ratios 

Starting from the weak amplitudes A^ defined in Section 2 and modifying them with 
FSI effects as explained in the previous Section, we evaluated the rates for all Cabibbo 
allowed two-body decays and for Cabibbo forbidden D^ and Df decays as functions of 
the parameters of the fit. These are ^, the parameters Wpp and Wpv of the annihilation 
contributions, the axial charges acu = ctcd and ttcs, the mass and width of the scalar 
Kq resonance. We use the values (m^, "i^, mg, rric) = (4.5, 7.4, 150, 1500 MeV) for 
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the quark masses. The decay constants are (/^r, fx, fp = fK*i fuji fcj)) = (133, 160, 
216, 156, 233 MeV). The pole masses in the form factors (|2.10 ) corresponding to yet 



undetected charmed particles have been taken to be M£)*(o+) = 2470 MeV and M£)«(o+) = 
2600 MeV. Other parameters relevant for decay to final states containing ri{ri') have been 
fixed following [^. The results are presented in Tables 1, 2, 3. 

The best-fit results are reported in column three of the Tables. The total y^ is 80.6. 
The 25 data points for Cabibbo allowed decays contribute 61.8 and the 12 data points for 
Cabibbo first forbidden decays 18.8. The best fit parameters are ^ = —0.027, r = —0.84, 
Wpp = -0.29, Wpv = +0.29, a^u = a^d = 1-0, a^s = 0.59, Mk* = 1928 MeV and 
Fk* = 300 MeV. A separate fit to Cabibbo allowed data alone gives quite similar values 
for the parameters. In the Tables we have also reported the theoretical predictions for 
the decays to final states containing K^, in order to show the importance of interference 
effects with doubly Cabibbo-suppressed amplitudes. 

We note that the relatively large SU(3) violation present in the data for exotic D^ 
decays [^ 



R 



+ 



Ur. 



Ucd 



2 



r(D 



+ ^ tt+ttO 



3.55, (4.1) 



T{D+ -^ KsTT+ 

(instead of i?_|_ = 1) is well reproduced by the fitted data, as it also happened in ^j. This 
point has been discussed in detail in [^ and more recently in [^. The reason is that 



several SU(3) breaking effects, each one rather small, contribute coherently to enhance 
R+. 

We have also evaluated, not including them in the fit, two recently measured doubly 
Cabibbo forbidden branching ratios. The experimental data are [^, [^ 

f.^ B.R.iD^ -^ TX-K+) , , . , ^ 

B.R.{D+^<t>K+) ,..+3.2 I „.^.„-2 ,4^^ 

^+- B.R.{D+^<Pn+) -M_2.6±0-7)10 ' (^.3) 

and the theoretical predictions using the best-fit parameters are i?^^ = 0.89 10~^and 
R'^ = 0.76 10-2. 

In the SU(3) limit the U-spin properties of the hamiltonian should give Rq^ = 
tan^ 6'c' = 0.26 10"^, also discussed in |J2^. Our model predicts Rq'^ to be 3.4 times larger 
than this value. This is due mainly to the W-exchange contributions, that should vanish in 
the symmetric limit and have opposite signs in the two amplitudes, and also to rescattering 
effects. 
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On the other hand the theoretical value for i?^ is much smaller than the experimen- 
tal datum, that however differs from zero by only 2.5 standard deviations. We note that 
in a factorized model the decay D^ -^ (pK~^ may only proceed through annihilation and 
rescattering. It is therefore difficult to reproduce the present very large value for R'_^_. 

We consider now D^ decay processes, and in particular Cabibbo first-forbidden 
decays. The D'^ meson is a t/-spin singlet and it should only decay to t/-spin triplet states, 
if fiavour SU(3) is a good symmetry. Therefore in that limit several relations among decay 
amplitudes hold. For the parity violating D^ -^ PP decays they are: 



AiD"" -^ K^'K ) = , (4.4) 

A{D^ -^ K+K-) = - A{D^ -^ TT+TT-) = 
^^(I>°^A-^0_^^^ (4.5) 

tan^c 
(the last equality corresponds to the symmetry prediction discussed after {K 



-^A{D^ -^ Ksn^) = cos e^^>A{D^ -^ Ksv) + sin e^^>A{D^ -^ Ksv') , (4.6) 
v3 

— ^A{D^ -^ 7r°7r°) = cos9^^>A{D° -^ w^rj) + sin 9 r^^>A{D^ -^ iv^r]') , (4.7) 

v3 

~A{D^ -^ 7r°7r°) = (1 - tan^ e^r,')A{D^ ^ r]r]) + 2 tane^r,'A{D^ -^ r]7]') , (4.8) 

A{D^ -^ rjTj') =tan^^^/ A{D^ -^ rjr]) + 

+ V3 {sme,^rj'A{D^ -^ n^T]) - cos e^^>A{D^ -^ n^r]')) . 

Comparing the above formulae to the experimental data 0,05 HH , we note that (|4.4|) is 
definitely not true, the moduli of the amplitudes in ( [4.5|) are in the ratios (1:0.59:1.15:0.67) 
instead of being equal, relation (|4.6|) is compatible with the data, but only with large 
phases, and finally no data exist for (^77|), ( [4.8|) , ( [4.9|) . 



(4.9) 



In the factorization scheme the amplitude A{D^ — > K'^K ) vanishes, and we may 
only obtain a nonzero rate in that channel through rescattering from the other decay chan- 
nels. Analogously, the smaU SU(3) breaking effects (^ 7^ /k, M|, - M^ > M|, - M|^,. . .) 
are not enough to reproduce the large ratio of K^K~ to tt+tt" rates. In Table 4 we show 
the results of our best fit to the Cabibbo first-forbidden decay rates. As explained in Section 
3, to the parameters previously determined we added two more free parameters, A and (f), 
and chose the value of another one, c, between the two solutions of a quadratic consistency 
condition. The fit we obtain is good (x^ = 1.7 with 4 data points). The best fit results for 

13 



the parameters are A = jMj, - Mj^ = 1205 MeV, (p = 48.7° and c ^ -2.69. The corre- 
sponding masses and widths of the scalar and isoscalar resonances are Mf^ = 1778 MeV, 
T/o = 361 MeV, M// = 2148 MeV and T// = 389 MeV. 

For parity conserving decays, the relations obtained from flavour symmetry are less 
restrictive in view of the few existing data. For instance, the relation corresponding to 
( [1.4| ), namely 

AiD^ -^ K^K*^) + AiD^ -^ K^K*^) = , (4.10) 



is valid in the factorization approximation and rescattering does not spoil its validity in 
our model. Only upper limits exist experimentally. 
Similarly, other SU(3) predictions are 

^(D° ^ K+K*-) = - A(D^ -^ n+p-) 

(4.11) 
A{D^ ^ K- K*+) = - A{D^ ^ n- p+) . 

The W-exchange terms strongly violate these relations and our predictions are therefore 
at variance with them. Data for np final states are unfortunately still missing. 

The experimental data and the fitted values for parity conserving Cabibbo first- 
forbidden D^ —>■ PV decays are also reported in Table 4. The only parameter added 
to those determined in the previous fits is the phase-shift of the isosinglet octet part of 
the decay amplitudes. This parameter turns out to be 243° and the quality of the fit is 
reasonably good (x^ — 7.7 for 4 data points). 

5. CP violation 

It is well known that CP violating effects show up in a decay process only if the 
decay amplitude is the sum of two different parts, whose phases are made of a weak (CKM) 
and a strong (final state interaction) contribution. The weak contributions to the phases 
change sign when going to the CP-conjugate process, while the strong ones do not. Let us 
denote a generic decay amplitude of this type by 

A = A e'^' + B e'^^ (5.1) 

and the corresponding CP conjugate amplitude by 

A = A*e'^'+B*e'^^. (5.2) 

The CP violating asymmetry in the decay rates will be therefore 

^ \A\'-\A\' ^ 2 -sjAB*) sini52 - 5^) 

^^"|^|2 + |^|2 \A\^ + \B\2 + 2^{AB*)cos{S2-5i)' ^ ' 
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Both factors in the numerator of eq. ( |5.3| ) should be nonvanishing to have a nonzero effect. 
Moreover, to have a sizeable asymmetry the moduli of the two amplitudes A and B should 
not differ too much. 

In Cabibbo first-forbidden D decays, the penguin terms in the effective hamiltonian 
( pr^ ) provide the different phases of the weak amplitudes A and B. Having obtained a 
reasonable description of the decay processes, including a model for their strong phases, we 
may envisage the more ambitious goal to derive CP violating asymmetries using our model 
for the phase-shifts. The asymmetries resulting are around 10"'^, somewhat larger than 
previously expected. We stress however that the actual numbers may vary appreciably for 
parameter variations that still give reasonable fits to the decay rates. 

For D~^ (and Df) decays the total charge allows to directly measure the rates to be 
combined in the asymmetry (|5.3|) . In the neutral D decays, however, the need of tagging 
the decaying particle to tell its charm and the possibility oi D — D mixing make (|5.3|) , as 
it stands, not directly measurable. Let us consider the particular case of an e+e~ factory 
at the -0", assume to tag D^ by semileptonic decays and define 



2//. _ l^(/»,e" 



\A{he 



+ ^|2 



"CP 



|^(/„e-)|2-M^(/„e 



+ 1|2 



(5.4) 



Limiting our considerations to the simplest case, D decays to a CP self-conjugate final 
state {fi = ±/i), and considering the time-integrated asymmetry, one obtains 



eff. 
'■CP 



1- 



1 + 



qA 
pA 



+ 



1 



y 



l + x' 



1 + 



1- 



qA 
pA 



1 + 



1 + 



qA 
pA 



+ 



l + x' 



qA 
pA 



(5.5) 



In ( |5.5|) A = A{D'^ ^ fi) is the decay amplitude, the mass matrix eigenstates are defined 



as 



\Ds,l) oc \D') ± ^ \d') 
p 



and the mixing parameters are 



X 



2 {Ms - Ml) 
Ts + Tr ' 



Ts + F, 



The mixing for charmed mesons is experimentally known to be small (|x| < 0.083, 
y < 0.085) and the theoretical calculations of the short-distance contributions give very 
small predictions. A reliable calculation of long-distance terms is problematic, even more 
so for CP violation in the mass matrix and the ratio p/q. We did not consider at all the 
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time-dependence in the asymmetries, since these depend on the phase of p/q. We note 
anyhow that the smaUness of x and y prevents the development in time of appreciable 
asymmetries even if the phases would allow this. We expect moreover that the modulus 
\p/q\ will differ from one by a small amount, O(10~"^). Therefore we can expand ( |5.5|) in 
the small quantities 



acp = 2 ( 1 



A 
A 



and 3fJe = - 
4 



1 



with the result 

a^^/- ~ acp (1 - 5) - 2 5 3?e + . . . (5.6) 

In (|5.6| ) we have used the definition 5 = 1 — (1 — y^) / (1 + x^). Experimentally 
5 < 0.014. Therefore, if acp is not much smaller than 3f?e, the measured asymmetry will 
be ttfjp ~ acp- 

We report in Table 5 the values of acp that we obtain in our model for several 
decay processes. We chose to give only the results that correspond to a good fit for the 
branching ratios, even if the predictions for some other decay channels are also large. For 
parity conserving D^ decays the final states are not CP eigenstates, and the corresponding 
formula for a^p is more complicated; however, we note that for amplitudes of not too 
different absolute values, as it happens in Cabibbo first-forbidden decays, and given the 
smallness of the mixing parameters, the result is again a^p ^ acp- If the final state 
contains a neutral kaon, normally identified by its decay to tt+tt", one has to disentangle 
the CP violating effects in D and K decays. How to do this for the D~^ decays has been 
discussed in [0. We note that these channels are not very promising candidates to look 
for CP violating effects in D decays. 

We evaluated the central value for acp choosing for the Maiani - Wolfenstein param- 
eters (p,77) the values (0.2, 0.3), following a recent analysis of CKM parameters [^ , and Ucb 



0.040. We varied [p,!]) in the one-sigma region obtained in [25] for fs = 200 ± 40 MeV. 



The error given in Table 5 refiects only this uncertainty, which is already quite large. 



6. Conclusion 

We have presented a generally successfull description of the complex of two-body 
nonleptonic decays of charmed mesons. We fitted 45 experimental branching ratios with 
eleven free parameters, that assume values close to the expected ones at the minimum of 
the x^ = 90. 
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We note that the large SU(3) breaking effects shown by the data are weU reproduced 
in our results. Rescattering (FSI) effects are particularly important in this respect: their 
parameters are derived from experimental data on nearby resonances. The masses of these 
being unequal, flavour SU(3) breaking is induced through the difference in the phase shifts 
for each isospin channel. 

W-exchange/annihilation contributions are substantial in many cases. The danger 
of getting too big decay rates for Dg Cabibbo-favoured decays has been avoided in our 
model imposing chiral symmetry requirements. 

Moreover, the rather large flnal state phase shifts and "penguin" operator contri- 
butions lead us to envisage CP violating asymmetries larger than 10~^ in some decay 
channels. 
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/. 


B.R.e.p(I50 ^ /.) 


B.R.tH{D^ -^ h) 


K-'K+ 


4.01 ±0.14 


4.03 


Ksn^ 


1.02 ±0.13 


0.78 


Kltt^ 


— 


0.57 


KsV 


0.34 ±0.06 


0.46 


KlV 


— 


0.34 


Ksi 


0.83 ±0.15 


0.84 


Klv' 


— 


0.67 


K*\' 


3.0 ±0.4 


3.49 


Ksp' 


0.55 ±0.09 


0.49 


Klp' 


— 


0.39 


K*-n+ 


4.9 ±0.6 


4.69 


K-p+ 


10.4 ±1.3 


11.19 


r% 


1.9 ±0.5 


0.51 


r%' 


< 0.11 


0.005 


Ks uj 


1.0 ±0.2 


1.12 


Klu 


— 


1.04 


Ks<P 


0.415 ±0.060 


0.42 


Kl4> 


— 


0.48 



TABLE 1 

Decay branching ratios in percent for Cabibbo allowed two-body D^ nonleptonic decays. 
In the first column the experimental data are reported (upper bounds are 90% c.l.) 0, in 
the second column the theoretical values obtained in the best fit. 
See text for further explanations. 
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/. 


B.R.,,p{D+ ^ /,) 


B.R.tH{D+ ^ /,) 


Ksn+ 


1.37 ±0.15 


1.08 


Kl7V+ 


— 


1.43 


rv+ 


2.2 ±0.4 


0.64 


Ksp+ 


3.30 ±1.25 


5.28 


KlP+ 


— 


6.49 


tt+ttO 


0.25 ±0.07 


0.17 


n~^r] 


0.75 ±0.25 


0.36 


Tx^rj' 


<0.9 


0.79 


k'k+ 


0.78 ±0.17 


0.86 


p07r+ 


<0.14 


0.17 


P+tt" 


— 


0.37 


p+r] 


< 1.2 


0.0002 


P+V' 


< 1.5 


0.13 


UJTl^ 


<0.7 


0.035 


(pTl^ 


0.67 ±0.08 


0.59 


k''k*+ 


— 


1.70 


k*'k+ 


0.51 ±0.10 


0.25 



TABLE 2 

Decay branching ratios in percent for Cabibbo allowed and first-forbidden two-body D~^ 

nonleptonic decays. 

In the first column the experimental data are reported (upper bounds are 90% c.l.)0, in 

the second column the theoretical rates obtained in the best fit. 

See text for further explanations. 
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/. 


B.R.e,^(D+ ^ /,) 


B.R.tH{Dt ^ /.) 


KsK+ 


1.75 ±0.35 


2.53 


KlK+ 


— 


2.26 


Tl^Tj 


1.90 ±0.40 


1.33 


n+Tj' 


4.7 ±1.4 


5.89 


p+T] 


10.0 ±2.2 


9.49 


p+v' 


12.0 ±3.0 


2.61 


k*'k+ 


3.3 ±0.5 


3.86 


KsK*+ 


2.1 ±0.5 


1.44 


KlK*+ 


— 


1.93 


(pn~^ 


3.5 ±0.4 


2.89 


uJ7r~^ 


< 1.7 


0.0 


P^-K + 


<0.28 


0.080 


p+n^ 


— 


0.080 


i^+TfO 


— 


0.16 


K+7] 


— 


0.27 


K+7]' 


— 


0.52 


K^n+ 


<0.7 


0.43 


K*+7rO 


— 


0.029 


K+p^ 


— 


0.24 


K*+r] 


— 


0.024 


K*+r]' 


— 


0.024 


K+u 


— 


0.072 


K+(j) 


<0.25 


0.015 


K*^7t + 


— 


0.33 


K^p+ 


— 


1.95 



TABLE 3 
Same as Table 2 for Df nonleptonic decays. 
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/. 


B.R.e,p(DO ^ /,) 


B.R.th{D^ ^ /,) 


n^T] 


— 


0.058 


7rO?7' 


— 


0.17 


rjr] 


— 


0.10 


rjrj' 


— 


0.22 


ttOttO 


0.088 ±0.023 


0.116 


TT+TT" 


0.159 ±0.012 


0.159 


K+K- 


0.454 ±0.029 


0.456 


k^k"" 


0.11 ±0.04 


0.093 


WTfO 


— 


0.008 


P'V 


— 


0.024 


pS 


— 


0.010 


UJl] 


— 


0.19 


tori' 


— 


0.0001 


(Pn^ 


— 


0.11 


<Pv 


— 


0.057 


j^*o^o 


< 0.08 


0.099 


rvo 


<0.15 


0.099 


K*+K- 


0.34 ±0.08 


0.45 


K*-K+ 


0.18 ±0.10 


0.28 


p'^Tl~ 


— 


0.82 


P~TX^ 


— 


0.65 


pOyrO 


— 


0.17 



TABLE 4 

Decay branching ratios in percent for Cabibbo first-forbidden two-body D^ nonleptonic 
decays. See the caption of Table 2 and text for details. 
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decay channel 


10"^ X acp 


decay channel 


10^ X acp 


D+ 


^p07r+ 


-1.17 ±0.68 


D+- 


.k'k+ 


-0.51 ±0.30 


D+ 


-^P+ttO 


+1.28 ±0.74 


D^ 


-^ttOt] 


±1.43 ±0.83 


D^- 


. K*0K° 


±0.67 ±0.39 


D^- 


-^ttOt]' 


-0.98 ±0.57 


D^- 


. k*'k^ 


±0.67 ±0.39 


D^ 


^r]r] 


±0.50 ±0.29 


D^^ 


K*+K- 


-0.038 ±0.022 


D^ 


^r]T]' 


±0.28 ±0.16 


D^^ 


K*-K+ 


-0.16 ±0.09 


D^- 


^^0^0 


-0.54 ±0.31 


1)0- 


-^ p^n~ 


-0.37 ±0.22 


D^- 


^ tt+tt" 


±0.02 ±0.01 


D^- 


-^ p~n^ 


±0.36 ±0.21 


D^^ 


K+K- 


±0.13 ±0.8 






D^- 


^k^k"" 


-0.28 ±0.16 



TABLE 5 

CP-violating decay asymmetries for some D^ and D^ Cabibbo forbidden decays. See text 
for explanation. 
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